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Abstract 

This paper concerns the diffusion-homogenization of transport equations when 
both the adimensionalized scale of the heterogeneities a and the adimensionalized 
mean- free path e converge to 0. When a = e, it is well known that the heteroge- 
neous transport solution converges to a homogenized diffusion solution. We are 
interested here in the situation where < e <C a <C 1 and in the respective rates 
of convergences to the homogenized limit and to the diffusive limit. Our main 
result is an approximation to the transport solution with an error term that is 
negligible compared to the maximum of a and ^. After establishing the diffusion- 
homogenization limit to the transport solution, we show that the corrector is 
dominated by an error to homogenization when a 2 < £ and by an an error to 
diffusion when e <C a 2 . 

Our regime of interest involves singular perturbations in the small parameter 
V = a' Disconnected local equilibria at rj = need to be reconnected to provide 
a global equilibrium on the cell of periodicity when r/ > 0. This reconnection 
between local and global equilibria is shown to hold when sufficient no- drift con- 
ditions are satisfied. The Hilbert expansion methodology followed in this paper 
builds on corrector theories for the result developed in [9]. 

1 Setting of the problem 

This paper studies the interaction between the convergence to a homogenized limit and 
the convergence to a diffusive limit in the context of linear transport equations, or 
linear Boltzmann equations [HIE], that model the propagation of particles in scattering 
environments. These two phenomena have been widely studied in the past [HI EH HI 
1231 I2T1 1221 1241 ITS] . The homogenization limit typically arises when the underlying 
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coefficients oscillate at a scale a <C 1 much smaller than the macroscopic scale at which 
phenomena are observed. The diffusion limit arises in highly scattering environments in 
which an equilibrium emerges in the velocity variable and a diffusion equation models 
the spatial behavior of the probability density of particles. High scattering is modeled 
by a mean free path £<1, where e is the (adimensionalized) main distance between 
successive collisions with the underlying medium. 

We assume here that the oscillations in the medium are periodic with period a < 1. 
When both a and e converge to 0, we expect the transport solution to converge to the 
solution of a homogenized diffusion equation. How fast convergence may occur and what 
are the main contributions to the error between the heterogeneous transport solution 
and its homogenized limit are the main problems of interest in this paper. We restrict 
ourselves to the case where the mean free path is (much) smaller than the correlation 
length a. We thus set rj = - and we assume that i)<l. The equation for the particle 
density takes the form: 

v ■ V£ + efl + -Q a fl = eS a (x, v) : = eS(x, -, v), (1) 
£ a 

in an infinite domain IR d x V where V is a smooth, symetric, compact domain in the 
velocity space such that ^ V. This equation could also be seen as an evolution 
equation by a change of variables u £ a = e~*/a in which the source S would play the role 
of an initial condition. We restrict ourselves to the time independent setting and to the 
case of constant absorption to simplify. 

The source term S is 1— periodic in y = -. We define the collision operator as 

Q a f = £(-, v)f(x, v) - [ a (-, v', v)f(x, v')dv{v'). (2) 
a J v a 

The choice of S is such that the above operator is conservative in the following sense. 
Define 

Qf(y, v) = E(y, v)f(y, v) - / a(y, v', v)f(y, v')du(v'), y G Y, (3) 

Jv 

where Y = [0, l] d is the unit cell. Then we assume the existence of C _1 > if) v (y,v) > 
C > such that 

(rjv ■ V y + Q)ip = 0, Y x V. (4) 

We also define 

{-r}v ■ V y + Q*)^* =0, Y x V, (5) 
where Q* is the adjoint operator to Q defined for a.e. y G Y and g G L 2 (V) as 

Q*g(v) = Y,(y,v)g(v) - / a(y,v,v')g(v')dv(v'). (6) 

Jv 

The simultaneous limit when the mean free path and the correlation length go to 
zero together with a = e has been considered in several recent papers [21 El IH El HOI El 
EH1 EH [201 [22] . In all those paper, only the case r\ = 1 is considered except in [23], where 
both ip v and tp v * are space independent and hence also independent of rj. A first result 
in the case r] 1 was obtained in [9] for spatially dependent ip v in the setting where 
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ijj*' 11 = 1. In that paper, a two scale convergence result is established in the setting 
where the limiting behavior involves the homogenization of a heterogeneous diffusion. 
The proof of that result is based on the weak formulation of the transport equation and 
on the method of moments. 

The aim of the present work is to implement the Hilbert expansion method, which 
leads to a strong convergence result under appropriate smoothness conditions in the 
(restricted) setting where if)* ,r> = i/j*(v). This strong convergence has two advantages. 
First, we are able to compute the correction terms in the expansion of the transport 
solution up to a term that is negligible compared to the maximum of a and 77 = -. 
Second, we expect this result also to be the first step in the study of the nonlinear 
Boltzmann equation given by the Fermi-Dirac equation in the spirit of [HE]. Note that 
in the case where 77 goes to infinity, we expect the limit to be the diffusion approximation 
to the homogenized transport equation obtained in [T5] . 

The main mathematical difficulty of the present work arises because the limit 77—7-0 
in (jlj) is singular. In the limit 77 = 0, local equilibria are obtained as a function of v £ V 
for each position y £ Y. For 77 > 0, a global equilibrium on Y x V emerges, as in the 
standard procedure obtained when 77 = 1 [21 HSl EEH EH E2J EE]- The passage from the 
local equilibria to the global equilibrium may in fact be arbitrarily complicated. Several 
conditions need to be imposed in order for a well-defined macroscopic equilibrium to 
arise. Even in the case 77 = 1 do we need to impose a no-drift condition. In the presence 
of drift, advection dominates scattering and entirely different phenomena arise (see, 
e.g., [32], [S] and the derivation of Euler equations when advection is dominant). Under 
appropriate sufficient symmetry assumptions similar to (though more constraining than) 
those in [2], we are able to verify the necessary no-drift conditions used in our derivation. 

The rest of the paper is structured as follows. The main hypotheses of regularity and 
no-drift as well as the main results of this paper are described in section El The main 
result is Theorem 12.11 below. The rest of the paper is devoted to its proof. Global a 
priori estimates are formulated and proved in section [3) The expansion of the transport 
solution in powers of s is treated in section HI The expansion in 77 of several cell solutions 
is given in section [5j The definition and 77-dependence of spatial density terms are given 
in section These results are combined to finish the proof of the theorem in section [7J 
Several technical results obtained in [5] and the generalization of their proofs if necessary 
are collected in the Appendix. 

2 Main result 

Under regularity assumptions recalled below, a standard application of a Banach fixed 
point theorem [13] ensures that (JTJ) admits a unique solution in H k (M. d , L 2 (V)) . The 
limit of f £,v (x,v) := f^(x,v) as 77 — > 0, however, involves singular perturbations. The 
reason is that the local equilibria in (J4j) and dSD become degenerate in the limit 77—7-0. 
In this limit, equilibria at different points y £ Y become disconnected and this can 
result in a very large effect at the macroscopic scale x. We consider here situations 
where the equilibria remain smooth in the y variable and generate no drift. Drift effects 
are ubiquitous in the homogenization of transport equations, with drastic effects as may 
be seen in, e.g., [5]. Diffusion limits arise under sufficient no-drift conditions as in, 
e.g., [2], which ensure that transport is not in an advection-dominated regime. Our 
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analysis in this paper shows that diffusion-like regimes are still valid in the singular 
limit ?7 — under appropriate no-drift assumptions. We show that these assumptions 
are consequences of symmetries of the transport coefficients, which we now define. 

Our first main assumptions is that ip* ,v (y, v) is independent of y and hence of rj. Note 
that Q in (j3J) is independent of rj and it is therefore not clear why non-trivial solutions 
would exist for all values of rj. Here, we assume that a.e. y G Y, we have an equilibrium 
xjj*(v) solution of 

Q*( 1 p*) = inV. (7) 

The equilibrium solution is allowed to depend on v but is independent of y G Y. When 
a(y,v',v) and H(y,v) are continuous functions bounded above and below by positive 
constants, then ip* can be chosen positive and normalized so that f v (ij}*) 2 (y, v)du(v) = 1 
a.e. y G Y. Moreover, up to normalization, ip*{v) is the unique solution to (j7j) as an 
application of the Krein Rutman theorem [21 [13] for the compact operators defined for 
a.e. y &Y: 



v7 r / a(y,v,v')f(v')du(v') 



v 



which preserve the (solid) cone of positive continuous functions. 

The no-drift conditions mentioned above will be verified under sufficient symmetry 
assumptions. We first assume that: 

a(y,v',v) = a(y,-v',-v) = a(-y,v',v), £(?/, v) = E(y, -v) = v). (9) 

We deduce that if)*(—v) = if)*(v). Note that a method to construct local equilibria 
consists of selecting a satisfying the above symmetries, if)* arbitrary (uniformly positive) 
such that v) = if)*(v), and finally define H(y,v) by ([7J), which also satisfies (jHJ). 

We also obtain the existence of a unique, bounded, positive, solution ijj(y,v) of the 
adjoint equation 

Q^(y,-)) = 0, inV, (10) 

normalized so that f v ip(y,v)ip*(y,v)dv(v) = 1 a.e. y G Y. It is not difficult to observe 
that ip(y, —v) = ip(y, v) = ip(—y, v) when (JSJ) holds. 

The Krein Rutman theorem for (JSJ) (all eigenvalues not equal to 1 have modulus 
strictly smaller than 1) shows that is a simple eigenvalue associated to Q and Q* and 
that all other eigenvalues of Q and Q* have strictly positive real part [21 [13]. On the 
vector space of functions / G L 2 (V) such that f v f(v)ip*(v)du(v) = 0, we define 



/"OO 

Q~ X f ■= - / e~ rQ fdr, (li; 
Jo 



which converges a.e. y G Y thanks to the spectral gap we just mentioned. Note that 
(Q -1 /, iP*)l 2 (v) = by construction. The inverse operator Q~* := (Q*)" 1 is defined 
similarly. We verify that under both Q~ x and Q~* preserve the subspaces of even 
and odd functions in the variable v. 

With tp* seen as a normalized solution of (— r/v ■ V y + Q*)ip* = on Y x V, we also 
obtain the existence of unique, bounded, positive, solutions ip v (y,v) of (JJJ) normalized 
such that fy xV if) v (y,v)if)*(v)dv(v) = 1. Upon introducing T v = rjv ■ V y + Q, we also 
define the inverse operator 

T^f = - f°° e~ rTV fdr, f G L\Y x V) s.t. / f(y, v)P(v)dydv{v) = 0. (12) 

JO JYxV 
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Finally, we observe that under ^(y^—v) and i/j v (—y,v) are also solutions of 
01]). Once properly normalized, since tp*(—v) = tf>*(v), we deduce that ^(y^v) = 
ip v (y, -v) = ip v (-y,v). 

Let us collect our main Assumptions: 

(HI) The velocity variables v lies in a compact, symmetric set V of M. d equiped with 
a symmetric probability measure v. There exist constants C, 7 > such that 
v({v e V, \v ■ £| < h}) < CW, for all f e S^" 1 , ft > . In particular, e 
V, |v ■ C| = 0}) = for all f ^ so that 

u • £ = a.e. in v implies £ = . 
(H2) The source term S G fl" 4 (R d (C°°(y), L 2 (V))) 

(H3) The scattering coefficient a is in C°(K x y„/;C~ r (Rj)). It is bounded from above 
and below by positive constants and is 1-periodic with respect to the variable y. 
The coefficient H(y,v) is defined implicitly in 

(H4) The symmetry relations OH]) hold so that the uniquely defined (after proper nor- 
malization) solutions ip*(v), i/j(y,v) and ij)' n {y,v) of (fTU|) . ((7|), and (@J, respec- 
tively, satisfy i/j*(—v) = ip*(v) and ip(y,v) = i/j(—y,v) = ip(y,—v) as well as 
^(y, = iF>(-y, v) = ip v (y, -v). 

We want to stress again that the limit 77 — ^ is singular. The limit of i/j 71 as rj — > is 
therefore not necessarily equal to tp(y,v). We are now ready to state our main results 
on the Hilbert expansion of the solution f £,ri := f£ of ([Q). The main result of this paper 
is the following theorem, which provides a strong convergence result for the corrector to 
homogenization theory: 

Theorem 2.1 Let f e ' r ' be the solution of (Tjp. Then the following expansion holds 



- v 



n 0,0 (x) g-l ( _ v . Vi/ ( p 0(JEty(!E „))) + n 0,l( x ) p ( ffi^ ( ffi j 



(13) 



77 L e J L 2 (R d xV) ' 

where ip(y,v) and ifi*(v) are solutions of ( flOl) and (J7j) ; respectively, and Q~ l is defined 
in ( 1TTT) . T/ie microscopic density p (y) is the unique solution of the elliptic equation 



L(p°) = £/ie normalization j p°(y)dy 

where 



Y 



L(P) = ~ / r(v)v ■Vy(Q- 1 (vV y (i;(y,v)p(y))))du(v). 
Jv 

The function 6~ x is given by 

9- 1 = L- 1 ( [ (vQ-\-v ■ V y p°(y)ij(y, v)) - v ■ V y (Q^(v^(y, v))p°(y)))^(v)du(v) 



v 



5 



with L 1 defined in Proposition \A.1\ below on functions in L 2 (Y) that average to on 
Y and the macroscopic density is given by the diffusion equation 

n°'°{x) - V x - (D-V*n°'°(x)) = [ [ S{x, y, v)if>*{v)dydv{v). 

Jy Jv 

The diffusion coefficient in the preceding equation is given by the expression 

D = / / (T\y,v)®vp\y)^(y,v) + 9*-\y)r(v)®vQ~\-vV y (p (y)ij(y,v)))dis(v)dy 
Jv Jy 

in which we have defined 

e*- 1 = L*~ x (j il>{y,v)v V y Q*-\vij*{v))dv(v)). 

where L* _1 is also defined in Proposition \A.l\ The correctors n 0,1 (x) and ri 1,_1 (x) satisfy 
the same elliptic equation as n°'°(x) with different source terms. Their expression is 
given explicitly in Proposition \6.2\ below. 

Before proving this theorem, we make a few remarks. 

Remark 2.2 The leading term in the expansion of f e a is given by n '°(x)p°(-)i(j(-,v). 
The behavior at the macroscopic level is given by n°'°(x), the solution of a standard 
diffusion equation. The microscopic level is given by the product of two terms. The first 
contribution to the product is the standard local solution ip(y,v), which indicates how 
particles are distributed in the v variable for each y G Y. The second, less standard, 
contribution is given by and indicates how the local (for each y) equilibria are 

related to one-another to generate a global equilibrium at the level of the cell Y. 

Remark 2.3 The above expansion implies that when 77 <C -, then the corrector is given 
by 

rV(y,«) • V,n°>°(x) + n^\x)p°(^)^,v). 

This is a regime of (relatively) low scattering where the corrector to homogenization 
(characterized by a term linear in Vn 0,0 ) dominates. The contribution n 1 " 1 provides a 
correction to the influence of the local equilibria at each y e Y to a global equilibrium 
on Y . 

When I 77, the corrector is given instead by 

n°'°(x)Q- 1 (-v V,(p°(^M^, «))) + n°< W(^M^, «). 

This is the regime of high scattering, where the correction to approximating the trans- 
port solution by a diffusion approximation dominates the correction coming from the 
homogenization procedure. The passage from local to global equilibria on Y generates 
a corrector described by n°' l (x). 

The rest of the paper is devoted to the proof of the theorem. 
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3 A priori estimates 

We start with an estimate that controls the remainder terms in the Hilbert expansion: 
Proposition 3.1 Let f £ a be the solution of (QJ). Then we have: 

WfJ + \\\fl - II < c(\\eS a \\ + ||S^f ||), (14) 

where || • || is the L 2 (M> d x V)-norm and for u G L 2 (R d x V) we have defined 

1 f 

u(x) = — — - / u(x, v)dv{v). (15) 



\V\j v 
Proof. We verify that 

(rjv V y + Q ri )l = 0, (-r]v V V + QT){W*) = 0, Y x V, (16) 

where we have defined the rescaled transport operator 

QMv,v)= [ g(y ' y. (y ' V>) [u(y, v) - u(y, v>)] dv(v>) . (17) 

Jv tp v {y,v) 

The reason for introducing the operator Q v is that 

(v ■ V x + -Ql + e)u = eF a {x, v) := eF(x, -,v), (18) 
e a 

where we have defined 

«) = 5fe4> v. «) = %^r> (is) 
ij r i{-,v) il>"(y,v) 

and 

Q> = / y ^(i^) [^«)-«(a,t/)]cfa/(t/). (20) 
Define the operator 

z2 = « • v, + (21) 

Then we recast (fTBT) as (e + T^)it = We verify that 
Ql(h) := (T2h,r a r a *h) 

= 1 / ^^^i^^ti^n^^) 1 ^'^"^ 1 '^ 12 ^)^)' (22) 



This comes from the fact that u • V^^*) = i<2a*(VW) so that 

e I 
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Let p = p(x) and s = u — p for an arbitrary p(x) independent of v. Then we 
find that Q^(p) = and more importantly that QP a {u) = Now for s such that 

J v sdu{y) = 0, we find that 

exoo > fikir, (23) 

for some f3 > as is clear from ( 122]) provided that tpa^T l § bounded from below by 
a positive constant uniformly in j]. We thus define p = u and s = -u — u so that ?2 is 
independent of v and s mean zero in v. Multiplying (jl8p by uip^T an< ^ integrating 
yields 

e(«, + = u ^X)- (24) 

As a consequence, 

f|M| 2 + ^IN| 2 < |(eF a , sr a r a *)\ + |(eF Q «>)|. (25) 
From this, we deduce the a priori estimate 

||u|| + -||«-«|| < c(||eF a || + \\F^X ||). (26) 
In the variables this is equivalent to (fT4|) . □ 



4 Expansion in e 

To emphasize the dependency in rj, let us denote f e,ri := f e a the solution of (CQ). As in 
the standard derivation of diffusion approximations, we first expand f e,v in powers of e 
at a fixed (arbitrary) value of rj. We prove the following proposition. 

Proposition 4.1 The solution f £,v can be expanded as follows 

f*n = n°> r '(xU''(— t v) + ef'^x, ^, v) + e 2 f 2 ' v (x, —,v) + £ 3 f' v (x, ^, v) + r £ '"(x, v) 
£ £ e e 

(27) 

where we have defined: 

f,v = T v-i(_ v . Vxn hv r)+ f* 



f 2,v = T v ~ l (S - v ■ y x T v ~ 1 {n°' r, il) v ) - J 0, 

The operator T ri ~ 1 is defined in ffl2|) . The density n°' v satisfies the diffusion equation 
n°^-W x - [ vx v {y>v)^*(v)diy(v)dy\7 x n^ = ! S{x,y,v)^*(v)du(v)dy (28) 

JYxV JYxV 

where x v — T* 7-1 (vtp v ) and the density for the correcting term n 1,r) satisfies 
n^-V x - fv X v {y,v)r{v)dv{v)dyV x n^ = - [v ■ V x f'\x, y, v)^{v)du{v)dy. (29) 

JYxV JYxV 



The remainder term r e,r? satisfies the following estimate: 



+s 2 ||/ 2 '1U^ i? + £ 3 ||/^||x^. 
Proof. We propose the following ansatz for / e '' 7 : 

We plug the ansatz into the transport equation and equate like powers of e to obtain 
the following sequence of equations. At the leading order, we obtain 

T r 'f°' r ' = which yields that f ' 71 = n '" (x)^ (y , v). 
We recall that T v = rjv ■ V y + Q. At the next order, f 1,v satisfies 

T v f i,v = _ v . V,/°'"(x, y, v) = -r(y, v)v ■ V x n^(x). 
We can then rewrite f 1,v as 

f^(x, y, v) = - X v (y, v) ■ V x n°"(x) + n^{x)r{y, v) 
by defining x v — T 11 ' 1 ^^ 11 ), the unique solution to 

T r> x r] = v ' l l )r> such that / / "x?(y, v)ifj*(v)dydu(v) — 0. 

Jv Jy 

We can apply the inverse operator T^' -1 defined in ()12p to the source term because 
the following standard no-drift condition is satisfied: 

v%l) r) {y,v) r il)*{v)du{y)dy = 0, 




v jy 



thanks to assumption (H4). 

The second-order equation reads 

T vf,v = _f,v {Xi v) + S ^ y iV )- v . VJ^( X , y, v) 

and the compatibility equation that the right-hand side must satisfy to be in the domain 
of definition of T v ~ l gives the diffusion equation (|28p for n 0,v . 

Introduce / ' the solution to 

jnff* = _ f o, v ^ ^ v) + ViV) _ v . v x (- x v (y, v) ■ V x n°«(x)). 
The third-order equation is the following 

TV 3 '" = -v ■ VJ^ix, y, v) - f^(x, y, v). 
and the corresponding compatibility equation gives a diffusion equation for n 1 '* 1 

n 1 ' v (x)-V x - vx v {y,v)ijj*(v)di>(v)dyV x n 1 ' r '(x) = - / v-V x f 2,r, (x,y,v)i>*(v)du(v)dy 

JYxV JYxV 
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since 

/ / r(v)x v (y,v)-V x n°'\x)=0. 
Jy Jv 

With the above expressions, we verify that remainder term now satisfies that: 

e e z ee e 

Thanks to the estimate (1141) proved in the previous section, we obtain (1301) . □ 

In order to analyze the behavior of the terms f ^, f 1,TI , f 2,r) and f 3 ' v as rj — > 0, we need 

to study some auxiliary equations. Such studies are conducted in the following section. 

5 Expansion of the auxiliary functions 

This section is devoted to the asymptotic expansion as 77 — > of the three cell functions 

i/j t >,X v , and x v *- 

5.1 Expansion of 1^ 

Recall that ip v satisfies 

rjv ■ Vytp v + Q{ip v ) = and / / ^{y, v)i/j*{v)du(v)dy = 1. 

Jy Jv 

We prove the following expansion 

Proposition 5.1 The function ifj v , solution to satisfies 

i] r >(y,v)=i) (y,v) + r]ip 1 (y,v) + r] 2 ij 2 (y,v) + rZ(y,v), H^H^^^) < Cn 3 , (31) 
with 

ip°(y,v) = p°(y)ip(y,v), 

^(y,v) = Q-\-v ■ V v ip°(y,v)) (32) 
tfj 2 {y,v) = Q~ 1 (-v-Vyij l {y,v)) + p 2 {y)^j{y), 

where 



L(p°) = with J p°(y)dy = 1 (33) 
p 2 = -L- 1 ( jf r(v)v ■ V y {Q~\-v ■ V y {Q-\-v ■ Vy^iy, v)))))du(v)^j . (34) 

The operator L and its inverse defined on functions in L 2 (Y) with vanishing average 
over Y are defined in Proposition \A.1\ in the appendix. From this expansion, we deduce 
that: 

W\\lIl~<C and || f vrrdv(v)\\ L ~<Cr). 

Jv 

Moreover, we find that -0° and ip 2 are even functions and if) 1 is an odd function in the 
variable v. 
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Proof. We plug the ansatz into the equation for ifj v and equate like powers of rj. The 
resulting series of equations is the following. First we have 

W°) = o, 

which implies that ip° = p°(y)if)(y,v). Next, we obtain 

Q(,l> 1 ) = -v-V v {p {y)il>{y,v)). 

Since thanks to assumption (H4) the no-drift condition f v vip(y,v)ip*(v)dv(v) = is 
satisfied, the compatibility condition for the previous equation is fulfilled and we can 
solve the equation for i/j 1 to obtain 

1> 1 = Q- 1 {-vV v {p°{y)*l>{y,v))). 

The compatibility condition for the third equation gives the equation for p° since for 

Q(^) = -vV y (^(y,v)) i.e., ^ = Q~\-v ■ V v ip\y,v)) + p 2 (y)^(y,v) 

to have a solution, we need 

/ iP*(v)vVy(ij 1 (y 1 v))dv(v) = 0. 
Jv 

This relation is satisfied if 

L(p°) = with L(p) = - [ r(v)v ■ V y (Q-\v ■ V„(V(j/, v)p{y))))dv{v). 

Jv 

The precise study of L and its adjoint L* defined by 

L*(p) = - [ ^J(y,v)v ■V y (Q*- 1 (vV y (r(v)p))du(v) 
Jv 

is conducted in the appendix. The next equation is 

The corresponding compatibility equation is satisfied since ip 2 is an even function of v. 
Then we write 

W 4 ) = -vV w (^ s ). 
Its compatibility condition gives the following equation for p 2 : 

HP 2 ) = ~ [ r(v)v ■ V y {Q-\-v ■ V y (Q~\-v ■ Vyifj^y, v)))))dv(v). (35) 
Jv 

Recall that L*(l) = 0. By the Fredholm alternative, equation (]35p thus admits a unique 
mean zero solution since 

/ / r(v)vV y (Q~ 1 (-vVy(Q- 1 (-vVyij\y,v)))))dv(v)dy = 0. 
Jy Jv 
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The remainder term satisfies 

rjv ■ Vr£ + Q{rf, = -rfv ■ V^ 4 ) 
and if we impose f p°(y)dy = 1 and J Y p 2 {y)dy = 0, we get 

/ / r%(y,v)if>*(v)du(v)dy = 0. 

JY JV 

Then using the technical regularity estimate (15 ip presented and proved in the appendix 
below, we obtain 

\\ r p\\Hk(Y),L^(V)) <V \\h1°+1(Y,L2(V))- 

The remainder involved in the proposition is given by r£ = r] 3 ip 3 + rfi\) A + and then 
satisfies 

ll^pll# fc (Y),L 2 (V0) < V 3 - 

By taking k large enough, we get the desired bounds. □ 

5.2 Expansion of x n 

We want to expand \ v satisfying 

rjv ■ V y x v + Q(X V ) = wJ> v and f [ x v {y, v)ifj*(v)dv(v)dy = 0. 

JY JV 

Proposition 5.2 We prove the following expansion 

X \y, v) = U-^yWy, v) + Q-\v^)p\y) - Q~\v ■ V,(r V(?/, v))) + r»(y, v), 
with | |f^| \H k (Y,L 2 {v)) — Cn, where we have 
L(p°) = with J p°(y)dy = 1 

LB' 1 = / (vQ-\-v ■ V yP °(y)ij(y, v))P(v) - wl>\v) ■ V y (Q-\v^)(y, v)p°(y)))du(v). 
Jv 

Moreover 

Wx n \\LlL T <- and || f v X n (y,v)r(v)dv(v)\\ L ~ <C. (36) 
V Jv 

Proof. In view of estimate (152]) obtained in the appendix below, we expect the 

expansion of x v to start at the order — . Thus, we begin with 

V 

Q( X ~ 1 ) = or x~\v,v)=6- 1 (y)il>(y,v), 
Then at the order 0(rj°) we get 

Q(X°) = W>° - v- Vyx- 1 - 
12 



Thanks to the no-drift condition and (H4), the compatibility condition for this equation 
is satisfied. Next, we write 

QOt 1 ) = W 1 - v ■ V yX ° i.e., x 1 = Q" W - v ■ V yX °)(y, v) + e^y^y, v). 

The compatibility condition for the equation giving y 1 gives the expression of 6" 1 since 
it satisfies the following elliptic equation 

L(0= / (vi; 1 (y,v)r(v)-vr(v)-V y (x(y,v)p (y)))dv(v). 
Jv 

This equation has a solution since 

/ / (vtp 1 {y,v)^*(v)-v^*{v)-V v (x{y,v)p°{y)))du{v)dy = 
Jy Jv 

because -0 1 defined in (j32j) is odd in y as a gradient of an even function and the second 
contribution can be written as a divergence in y and hence averages to over the cell 
Y . Note that 9~ l is odd in the y variable and that therefore 

X~ l (y,v)^*(v)dv(v)dy = 0. 

Finally, we write 

Q( X 2 ) =vj> 2 -v WyX 1 

for which the compatibility condition holds since if) 2 is even with respect to v and x 1 is 
even. Indeed, ip 1 is odd with repect to v (even with respect to (y,v)) and hence so is 
X°- 

The remainder term satisfies 

rjv ■ V„rjJ + Q(rJ) = -rfv ■ V yX 2 + vr%. 

We chose 9 l such that J v (—r/ 3 v ■ V y x 2 (y, v) + vr^(y,v))t()*(v)du(v) = for any y EY. 
Note that this is possible since J v f Y (~ r f v ' ^yX 2 + vr^)ip*diy(v)dy = because 
is even with respect to (y, v) whereas x 2 is odd. The fact that x 2 an d vr p are °dd 
with respect to (y, v) implies, thanks to Proposition IA.lt that 9 1 is odd and therefore 
Iy Iv X 1 (Vj v)i/j*(v)du(v)dy = 0, which finally gives 

/ / r^(y,v)ip*(v)du(v)dy = 0. 

By using estimate (15T|) and since || — • V^y 2 + w2||H fc (y;i 2 (v)) < C?? 2 ) we have that 

This concludes the proof of the result. □ 
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x*°(y,v) 


= Q* 


L*(9*- 1 ) 


-L 






L*(6*°) 


-L 






x*\y,v) 


= Q* 


L*{6* 1 ) 


-I 







5.3 Expansion of x 11 * 
Proposition 5.3 Let x v * be the solution of 

- V vv yX v * + Q*(x v *) = vr- 

The following expansion holds 

x v * = l -e*-\yW{v) + x *°(y, v) + vx*\y, v) + r^(y, v) 

with ||r^*||H fc fy,L 3 (v')) < Crf , where 

+ v . v^*-V*)(y, v) + e*\ y )r(v) (37) 

tMvVyQ'-^vipWduiv), (38) 

i;(y, v)v ■ V y Q*-\v ■ V.Q*" 1 ^* + v ■ V yX *- x ))dv{v), (39) 
^(v-VyX^y^ + e^r^v), (40) 
My, v)v ■ V y Q*-\v ■ V y Q*-\x*°))dv{v). (41) 

Proof. The proof is a straightforward adaptation of the proof of the same result in 

[9]. We briefly sketch it to be self-contained. As in the expansion of x v , we start at the 

order - and write 
v 

X 1 = -X +X +VX + V X + V X +r'* 

with 

Q*( X *^) = so that x*" 1 ^) = e*-\y)r{v). 

Then we have 

Q*(x*°) = vr + v-v yX *- 1 

which admits, thanks to the no-drift condition, a solution given by 

X *°(y, v) = Q*~\vr + v ■ V.x*" 1 )^, v) + 6*°(y)r(v). 
The first order equation is 

Q*(x* 1 )=v-V yX *°. 
The compatibility equation for this equation gives 

L*(8*~ 1 ) := — / My,v)v-V v Q*- 1 (v-V y (i(>*(v)0'- 1 )du(v)= [ ifj(y,v)vV y Q*~\v4j*(v))du(v) 
Jv Jv 

This equation has a solution since 

p\y) [ ifj(y,v)vV y Q*- 1 (vr(v))dv(v)dy = 

Y JV 
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because p°ifj is even with respect to the variable (y, v) whereas vip* is odd. Then x* 1 
can be written as 

x*\y, v) = Q*-\v ■ V yX * (y, «)) + e*\y)r(y, «)■ 

The compatibility condition for 

gives the equation determining 9*° since we need 

L*(6*°)= [ ^{y,v)vV y Q*-\v-V y Q*-\vr + v-V yX *- l ))dv{v). 
Jv 

Note that the parity properties ensure that the previous equation has a solution. The 
final equation is 

Q*(x* 3 ) = v-V yX * 2 
whose compatibility condition gives the equation for 9* 1 

L*(9* 1 ) = f ^(y,v)v -VyQ^iv -VyQ^ix^dHv). 
Jv 

The remainder term satisfies 

rjv ■ V w r£. - Q*{r v r ) = rfv ■ Vx* 3 - 

As for x* 7 , 9* 1 are odd for i = —1,0,1 and then, J y J Y r^*il}*dv{v)dy = 0, and then, 
thanks to estimate (155]) . we get 



\ r x*\\H k (Y,L 2 {V)) < Crj 1 



which concludes the proof. □ 



6 Estimates on the densities 

We are now ready to analyze the expansion in n of the densities n h,ri for k — 0, 1 obtained 
in Proposition 14. 1L We start with the following simple lemma: 

Lemma 6.1 Let n v be the solution of 

ra"(x) - V, • (D n ■ V x n v (x)) = S v (x) 

with 

jyq _ D° + rjD 1 + 7]D V , where \ \ \D V \ \ | — 0, || ■ ||| a norm on d x d matrices 
= S° + rjS 1 + V S V where S° G H k (R d ), S 1 G H k (R d ) and IMU^ -^o 0. 

We decompose 

n v = n° + rjn 1 + r)fi v 
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with 

n°(x) - V, • (D°V x n°(x)) = S°(x), \\n°\\ Hk+ i m < C\\S°\\ Hkm 

and 

n 1 (x)-V x iD°V x n 1 (x)) = S 1 (x)+V x -(D 1 -V x n°(x)), 1 l^l \ Hk+Hm < C\ \S°\ |^ (Rd) +| {S 1 
Then n v satisfies 

| \h v \ |_fffe+i(Rd) — > when 77—7-0. 

Proof. Indeed h v satisfies the equation 

n" - V, ■ (D°V x h ri ) - V, • [(77D 1 + r/D") V x n"] = 5" + V, ■ [D^V x n° + { V D l + • n D r >) V.n 1 ] . 

A standard a priori estimate shows that H^H^-i^d) — > when 77—7-0. Upon differenti- 
ating the equation k times, we obtain the result in H k+1 (WL d ). □ 

We apply this lemma to derive a convergence result for n ' v and n ' v . 

Proposition 6.2 Assume S G H k (R d ). First, we have that 

luo,77 ^0,0 „„o,iii . n 

I \n — n — 7777- 1 |H fc + i (R d ) ? u 

where n 0,0 G zs i/ie solution to 

n°'° — V x • (D • V^n ' ) = / / S(x,y,v)r(v)dydu(v) 

Jy Jv 

and n ' 1 G zs £/ie solution to 

n o,i — V x • (D • V^ ' 1 ) = V x ■ (D 1 • V x n°<°), 



y */y 



</y </y 

= / / <?~\v^* + v- V v F- 1 V){y,v)®vp\y)il>{y,v) 

JV JY 

+ fl*- 1 ^)^) ® vQ~\-v ■ V y (p°(y)i>{y, v)))dv{v)dy 

D 1 = / [ (X*~ 1 (y,v)®vip 2 (y,v) + x*°(y,v)®vip l (y,v) + x* 1 (y,v)(g)vip (y,v))dis(v)dy 
{6*-\y)r{v) ® vQ~ x {-v ■ V y {Q-\-v ■ V y (p\y)^y, v)))))dv{v)dy 
[Q*-\v^j* + v ■ VyO*- 1 ^*) + 0*V] ® • V y (p°(y)ij(y, v))du{v)dy 

Q*~\v ■ W v Q*-\vr + v ■ V y 0*- V) + 0*\y)r(y, «)) ® vp (y)ijj(y, v)du(v)dy. 

Wn 1 ' 11 - ^n 1 '^ 1 1 | H fc-i( R d) -^0 (42) 




y Jy 




y Jy 




'y Jy 

Second, we have that 



1(3 



with n 1 ' 1 G H k 1 (R d ) the solution of 

n 1 '- 1 - V x ■ (D ■ V^n 1 '- 1 ) = S ,1 '~ 1 (x) (43) 
where S ,1,_1 G iJ fc_2 (IR d ) groen by 

S^~\x) = -vJ [-n°'°(:z)^(^)+Stx,^)+^ 



JYxV 

-V x ■ [ i-v- V x (-x-\y,v) ■ V x n > (x))}x*°(y,v)dv(v)dy. 

JyxV 

Proof. Recall that n 0,ri satisfies 

n 0,r? + V x • (D v ■ V x n ' v ) = [ [ S(x,y,v)r(v)du(v)dy. 

Jy Jv 

By using the definition of D v = X v *{Vi v ) ® vip v (v)dydu(v) and the expansions of 

Jv Jy 

ip v and x v * obtained in the preceding section, we obtain that 

D v = D + 77D 1 + r]D v with HIZ^HI -»„_>o0. 
Then Lemma [6.11 gives the first result. On the another hand, n 1 ' 7 ? satisfies 

n 1 '" - V x . • (D v ■ V x n^) = S 1 '"^) 

with 

S^(x) = -V x - I I (f'\x,y,v))vr*(v)du(v)dy 
Jy Jv 

= -V, f [ T' V T«* X r >*(y,v)dv(v)dy = -V x [ [ T^f'\^(y,v)du(v)dy 
Jy Jv Jy Jv 

JYxV 

= -S^\x)+o(-). 
rj rj 

Since S 1,r * is of order -, we apply Lemma IBTTl to rfn l,v to obtain ( l42l) . This concludes 
the proof of the result. □ 



7 Proof of Theorem 2.1 



We are now in a position to conclude the proof of our main result. We first obtain an 
estimate for the remainder term r e,r? in (1271) . 

Proposition 7.1 The different terms involved in the estimate of the remainder term 
r e ' v are bounded as follows 

\\f°«\\ HiL 2 L ~<C (44) 

ll/ M ll^- W <^ and ||«.V,/^|| fl2 _ fcw <-g r l<fc<3. (45) 
These result yield that 
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Proof. Let us start with / 0,r? . Recall that f 0,7] (x,y,v) = n ' 71 (x)^ 71 (y , v) and that 

W\\lIl~<C and \\n^\\ m < C. 
This gives (PHI) . Concerning f l,v = —x^iViV) ' V x n 0,7] (x) + n 1,v (x)ifj ri (y,v), we get 



H/^llflSLSzs- < C[||n ^||^|| x "|| w + \\n^\\ H s\\r\\ L ^ y 

and 



\\v-V x f l ^^*\\ H 2 L 2 L ^ < \\n ^\ \ H 4 1 1 \vx v ip* 1 1 HlLlh^ + ||« ' v \\hi\\v^*\\hIlIl^ < C. 

This yields ( |45|) for = 1. Let us now consider / 2,r? . Since J y f Y f 2,T, ^)*{v)dv{v)dy = 0, 

IM/ 2,?? )lkw < | and || / v T"(/ 2 ^)V*(y)M^)^ll^ ^ C ' we § et the first P oint 
of P5|) for k = 2. Indeed, we decompose 



v 



and use (15T1) and (IB^i) . respectively, obtained in the appendix below. 
Concerning the second point, we write after defining x* — Q* -1 ^^*)-, 

v ■ VJ 2 '"r = / v-V x f 2 ^*diy(v) = [ V x f 2 ^-Q*(x*(y,v))du(v) 
Jv Jv 

= [ V x Q(f 2 «)- X *(y,v)dv(v) 
Jv 

= I V x (- V vV y f 2 ' r >-v-V x f 1 «(x,y,v)-f« + S)-x*(y,v)du(v), 
Jv 



so that 



x v y J 

< 



\\v ■ v x p^\\hilil T < c[ V \m H 2 L ^ + \m mL 2 LT + n/ o '"ii ff | W + \\s\\ mLU 

is bounded by Crf x . In a same way, since J v J Y f^^^du^dy = and | |T T? (/ 3,r ')| \hil%l 
|and || jyT^if^^^dv^dvWuiL^ < |, we get the first point of (0SD for k = 3. 
Finally, we write 

v ■ V x f 3 ^* = [ vV x f>^*dv{v) = [ V x f«.Q*( x *(y,v))dv(v) 
Jv Jv 

v x Q(f*)-x*(y,v)dv(v) 

V x {-7)v ■ V y f« - v ■ V x f«{x, y, v) - Z 1 '") ■ X *(y, v)du(v), 



so that 



\V V X P^*\\LILIL T < C[ V \\f«\\ HiLU ~ + ||/ 2 '1|^ W + H/^lkw 



c 

< — 

rf 
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which is bounded by Crj 2 . This concludes the proof of (|45p for k = 3. 
In view of the estimate (130]) and (1451) . we get 

This concludes the proof of the proposition. □ 

Now, the rest of the proof of the main theorem simply consists of collecting the 
expression for f £,v given in ( 127)) . the above estimate for r e,r? , the expansions for the den- 
sities in section and the expansions for the auxiliary functions ip v and x v — T^~ x (vi/j 71 ) 
given in section This yields (fT3")) and concludes the proof of Theorem 12.11 

8 Postscript 

Guillaume Bal and Marjolaine Puel concluded this work after Naoufel Ben Abdallah 
passed away the Fifth of July, 2010. We keep the memories of beautiful discussions we 
started in Montreal and continued in Toulouse. 

A Previous results 

In this appendix, we improve useful propositions proved in [9]. Their proofs in the case 
where ip* depends on v is a straightforward adaptation of the proofs obtained in [9J 
when ip* = 1, except for the estimates ()5Tj) and (152)) below which are new. 

Proposition A.l Let L be the unbounded operator on L 2 er (M.y) with domain H 2 er {W^), 
defined by 

L(p) = - [ r{v)v-V y (Q- 1 {v-V y {i(;(y,v)p(y))))du(v), 
Jv 

= -div y {DV yP ) + V y -{U (y)p) 
and let L* be its adjoint defined by 

L*(n) = — div y (D J (y)V y n) — U(y) ■ V y n with periodic boundary conditions. (46) 
The matrix-valued function D and the vector field U are given by 

D(y)= t il)*v®Q- l (vil>)dv{v), U(y)= / <g> Q~ l (v ■ V y ^)du(v). 
Jv Jv 

The following statements hold: 

(i) Im{L) = {ue L 2 per (R d y ), s.t. f Y udy = 0} . 

(ii) There exists a unique, nonnegative function p(y) 6 H 2 er {Ry) such that L(p) = and 
fypdy = 1 . 

(Hi) For any function v e Im(L), there exists a unique solution to L(u) = v, J Y udy = 
with u G Hp er (M.y). This solution will be denoted u = L~ l {v), and L~ l will be referred 
to as the pseudo inverse of L. 

(iv) Im(L*) = {ue L 2 per (R d y ), s.t. J Y u(y)p(y) dy = 0} . 
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(v) The kernel of L* is the set of constant functions (in the variable y). 

(vi) For any function v* £ Im(L*), there exists a unique solution to L*(u*) = v* , 
j Y u*(y)dy = with u* £ Hp er (Ry). This solution will be denoted u* = L*~ l (v*), and 
L*^ 1 will be referred to as the pseudo-inverse of L*. 

(vii) D(y), U(y) and p(y) have the following symmetry properties: 

• The diffusion matrix y > D(y) is even. The flux y i— > U(y) is odd. The equilibrium 
function y h-> p(y) is even. 

• The sets of even and odd functions in y are invariant under L _1 and L*~ x . 

(viii) The diffusion matrix D{y) is in C^ r (Ry) and satisfies 

D{y)£ ■ £ = f = and furthermore D(y)£ ■ £ > fi\£\ 2 , (47) 

for some constant (3 > 0. 

(ix) There exists a constant C such that 



|L*-V||_h-2 <C||u*|| L 2 , \\L- l v\\ H 2 <C||d| L 2 

\L*- l V*\\ L 2 <C\\v*\\ rT -2, ||L-M L2 <C\\v\\ f r- 



(48) 



2 



for all v £ Im(L) and v* £ Im(L*). 



Remark A. 2 We need the fact that ifj* does not depend on y to identify zero as the 
first eigenvalue of L* and then of L. When ip* is allowed to depend on y, more complex 
global equilibria (or possibly the lack of such a global equilibrium) need to be analyzed on 
Y. 



Proposition A. 3 The operator T v = wv ■ V y + Q[y](-) is an unbounded operator on 

2 (Tad 



Ll er (Ry x V) with domain 



V = {u£ L 2 per (R d y x V), such that v ■ W y u £ L 2 per (R* x V) } . 
Let T 11 * = -nv ■ V y + Q*[y}(-) be the adjoint ofT^. Then 

(i) The kernel of T v is a one dimensional space spanned by a positive function ip v (y, v) 
normalized by J VxY ilJ v ifj*(v)diy(v)dy = 1. The function ifj v satisfi 



es 



^(—y, —v) = ip v (y, v), a.e. in y and v . (49) 

(ii) The range of T v is the set of functions g £ Lp er (Ry x V) such that we have 
JvxY 9(V, v)ip*{v) du(v)dy = . 

(Hi) The adjoint T v * has the same domain T>. Its range is the set of functions g such 
that J VxY i> v (y, v) g(y,v)du(v)dy = 0. Its kernel is spanned by tp*. 

For g £ L 2 er (Jty x V) satisfying j VxY g(y, v)i/j*(v)dv(v)dy = 0, there exists a unique 
function EJ 1 £ T> such that 

T r >(K n )=g and [ R v (y,v)ip*(v)du(v)dy = 0. (50) 

JVxY 
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We denote RP = (T v ) 1 (g). There exists r] Q > such that 

C 

WR^Wh^Y^HV)) < ^\\g\\Hk(Y,L2(V)) , o<v<vo- (si) 
If in addition, J v gip*(v)dv(v) = 0, then 

C 

\\R v \\h><{y,l2(v)) < — \\g\\m{Y,LZ{v))- (52) 
Moreover, the following symmetry implications hold true 

Ifg(-y,-v)= g(y,v) a.e. then RP(-y,-v) = R v (y,v) a.e. 
If 9(-y,-v)= ~g(y,v) a.e. then RP(-y,-v) = -RP(y,v) a.e. 



(53) 



(iv) For g* G Lp er (Ry x V) satisfying f VxY ip v (y,v)g*(y,v)di>(v)dy = there exists a 
unique function RP* G T> such that 

T v *(R v *)=g* and ! R>*( y y,v)'ilj r >(y,v)dv(v)dy = . (54) 

JVxY 

We denote RP* = (T* 7 *)" 1 (g*). There exists r] > such that 

C 

\\R v *\\hx(y,l2(v)) < ^\\g*\\m(Y,tf{v)) , < rj < ?7o . (55) 
// moreover, j v g*ipdv(v) = 0, then 

C 

\\H k (Y,L2{V)) < — \\g*\\m{Y,L 2 {v))- (56) 
Finally, the following symmetry implications hold true 

Ifg*(-y,~v)= g*(y,v) a.e. then RP*(-y,-v) = R n *(y,v) a.e. 

ifg*(—y>— v )= —g*{y, v ) a - e - tnen R v *(—y,—v) = —R v *(y,v) a.e. 



(57) 



Proof. Step 1: For the proof of estimate ( 1521) . we argue by contradiction. Assume 
that 1 1^| | = 1 and ||^|| ->^ with 

nv ■ V y RP + Q(RP) = n~g\ 

Decompose RP = ^{y^ijj, v) + 5 v (y, v), we get by the dissipation property 

\\nh(Y,v) < C [ [ Q{fr)^d»(v)dy 



Y JV 



< [ I \vg v -Vv-V y R v )^*dv(v)dy 
Jy Jv y 

< C(77||^|U 2(yxy) ||^|U 2(yxy) +v\\R v \\h(YxV)) 
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which implies that | l^l \l 2 (yxV) ~~ when 5 goes to zero. 
On the other hand 

rjv ■ V y 5 r > + Q(cf') = r)g v - rjv ■ V^V) 
and an integration with respect to v against if)* gives 

divj, J vS^dviv) = J g^duiv) = 0. 
Recalling that x* = Q*' 1 ^^*), we write 

div y / v^dviv) = -div y ( [ X *Q(6 r >)dv(v)) 



Jv JV 

= rjdiv y ([ x*v-V v 8*dv(v))+7idiv v (f X *v-V y m))- V div y ([ X *9 v dv(v)) 
Jv Jv Jv 

which leads to the following elliptic equation for 7^ 

L( 7 ") = div„( / x*v-V y 5^du(v))-dw y (f X *~9 v dv(v)). 
Jv Jv 

Since the right-hand side goes to zero in H~ 2 (Y), we obtain that 7^ — C v p°(y) — > 0. 
But 

C- f l n dy = - [ L-\div y (f x*vV y 5Uv{v))-diw y {f X * ~9 v dv{v))) 0. 
Jy Jy Jv Jv 

Moreover 

l v (y)dy = [ <f{yMy,v)il>*{v)dydu{v)= [ {W{y) - 5 r '(y,v))r(v)dydu(v) 

JvxY JvxY 

8' n {y,v)ip*{v)dydv{v) — > when 77 — >■ 0. 

VxY 

Finally, we obtained that R v — > strongly which leads to a contradiction. 

Step 2: Proof of estimate (I5ip . We first recall that ||-R ,? ||l 2 (y,v) < C. A direct 
application of the estimate proved in [9] would make us lose two powers of r\ for each 
derivative. The proof of flBTl) requires additional computations. 

We prove the result by induction. Assume that for any multi-index i satisfying 
|i| < \k\, we have < C. To prove that \\dyR v \\ L 2^ Y y) < C, we argue by 

contradiction. 

First, we note that for any / and any multi-indice k = (kx, k^, • • • , kd), we have 

|fcl| \k d \ 

d k y (Q*(f)) = E • • • E nf=i(cg)^-^ • • • afc**Q{sfc ■ ■ ■ aft/). 



Thus, we write 



|fc|-i 

V v ■ V y d k y W + Q(tfRP) = v 2 d y g v - E QtfTQi&yW 

1*1=0 



22 



where = Ilf =1 (C^) and where the sum runs over i = (ji, 12, ■ ■ ■ ,id) such that |i| < 
\k\ — 1 and %i < k\ for 1 < I < d. To obtain a contradiction, we choose a renormalization 
that implies || c^i^H = 1 and ||9^' 7 || — > and ||£^J2 , '|| — > for any i satisfying |i| < 
We decompose dyffl into two parts as follows 

We obtain from the dissipation property that 



f f d k RP 

\\%\\h<Y,v) < ~C j Y J v Q(d k y W)^—rdv(v)dy 

= cf /(^X^-Ec^^ 

JyJv |i|=0 V 

which gives 

I \8k\ \l 2 (y, v) — V>o 0. 
On the other hand, 7^ satisfies the following elliptic equation 

Hit) = I dffir + div y ( [ X *v-VM))-Vdiv y [ affix* 
Jv Jv Jv 

+-div„ / Cid^QidlR^x*. 

We have to control the last term, which is a priori of order -. For that purpose, we 
write 

Jv Jv 
Moreover, since if)* does not depend on y 

Ifeil 

d k y (Q*(x*)) = E • " " E nf =1 (^;)djr- • ■ ■ sfc+Qift ■ ■ ■ d y *x) = 

tl id 

and then 

|fe|-i |fe|-i \k-i\ 

/ J2 cid^Qi&R^x* = - / e E c i-i d i^ j Q*( d jx*) 

Jv |i|=0 V i=0 \j\=l 

which after rearranging the sums gives 

r 1*1-1 f \k\ \k-j\ 

/ Cld^Q^R^x* = / E E cua k d k y ^Q(d l &) 

|*|=o b|=i |*|=o 

|fc| |fc-i| 



/ Eci^E^-j^ - ^^^)- 

" V 1*1 1 I * I n 



|il=i KI=o 
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Finally, this leads to 




This term converges to zero in H 2 (Y, L 2 (V)) when 77 — y and leads to a contradiction. 
□ 
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